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ABSTRACT 

Millimeter wave Very Long Baseline Interferometry (mm- VLBI) provides access to the emission region sur- 
rounding Sagittarius A*, the supermassive black hole at the center of the Milky Way, on sub-horizon scales. 
Recently, a closure phase of 0° ± 40° was reported on a triangle of Earth-sized baselines (SMT-CARMA-JCMT) 
representing a new constraint upon the structure and orientation of the emission region, independent from those 
provided by the previously measured 1.3mm- VLBI visibility amplitudes alone. Here, we compare this to the 
closure phases associated with a class of physically motivated, radiatively inefficient accretion flow models, and 
present predictions for future mm- VLBI experiments with the developing Event Horizon Telescope (EHT). We find 
that the accretion flow models are capable of producing a wide variety of closure phases on the SMT-CARMA- 
JCMT triangle, and thus not all models are consistent with the recent observations. However, those models that 
reproduce the 1.3mm- VLBI visibility amplitudes overwhelmingly have SMT-CARMA-JCMT closure phases be- 
tween ±30°, and are therefore broadly consistent with all current mm- VLBI observations. Improving station 
sensitivity by factors of a few, achievable by increases in bandwidth and phasing together multiple antennas at in- 
dividual sites, should result in physically relevant additional constraints upon the model parameters and eliminate 
the current 180° ambiguity on the source orientation. When additional stations are included, closure phases of 
order 45°-90° are typical. In all cases the EHT will be able to measure these with sufficient precision to produce 
dramatic improvements in the constraints upon the spin of Sgr A*. 

Subject headings: black hole physics — Galaxy: center — techniques: interferometric — submillimeter: general 
— accretion, accretion disks 
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1. INTRODUCTION 

It has recently become possible to study the emission regions 
of a handful of black holes on sub-horizon scales with millime- 
ter wave Very Long Baseline Interferometry (mm- VLBI). Al- 
ready, this technique has produced horizon-scale information 
on Sagittarius A* (Sgr A*), the 4.3 x 1O 6 M black hole located 
at the center of the Milky Way (Ghez et al. 2008; Gillessen et al. 
2009b, a), using the Arizona Radio Observatory Sub-Millimeter 
Telescope (SMT) on Mount Graham, Arizona, James Clerk 
Maxwell Telescope (JCMT) and Sub-Millimeter Array (SMA) 
atop Mauna Kea in Hawaii, and the Combined Array for Re- 
search in Millimeter-wave Astronomy (CARMA) in Cedar Flat, 
California (Doeleman et al. 2008; Fish et al. 201 1). Due to the 
limited signal-to-noise of these early experiments, they have 
produced primarily visibility amplitudes, related to the modu- 
lus of the Fourier transform of the intensity distribution of the 
source. Nevertheless, when analyzed in the context of phys- 
ically motivated accretion flow models, they have resulted in 
dramatic constraints upon the spin orientation and magnitude 
(Broderick et al. 2009; Huang et al. 2009; Dexter et al. 2010; 
Broderick et al. 2010). However, the absence of phase informa- 
tion introduces fundamental degeneracies in the orientation of 
the modeled image, and systematic uncertainties in the image 
structure generally. 

The importance of visibility phase information has been ap- 
preciated since the beginning of radio interferometry. In the 
context of VLBI, it has only been widely possible since the 
introduction of "closure phases", and the associated develop- 
ment of self-calibration techniques (e.g., Jennison 1958; Jenni- 
son & Latham 1959; Rogers et al. 1974; Pearson & Readhead 



1984; Thompson et al. 2001). The closure phases, which are 
equivalent to the argument of the bispectrum, are combinations 
of the visibility phases on baseline triangles (and thus VLBI 
station triples) that are insensitive to individual station-based 
phase errors, which otherwise typically dominate the phase un- 
certainties, e.g., due to pathlength variations from atmospheric 
turbulence. 1 Specifically, if 0, ; is the visibility phase on the 
baseline between stations ; and j, the closure phase associated 
with three stations is 



®ijk = 4>ii+4>jk + 4>M- 



(i) 



Even a handful of closure phases are diagnostic of the under- 
lying image structure, e.g., the closure phases of a point source 
or Gaussian flux distribution (including asymmetric cases) are 
identically 0°, while that of an annulus may be either 0° 
or 180°, depending upon the particular baselines considered. 
More general flux distributions produce non-vanishing closure 
phases, indicative of the symmetry of the image. Recently, Fish 
et al. (2011) reported the measurement of a closure phase at 
1.3 mm for Sgr A* of 0° ±40° on the SMT-JCMT-CARMA 
triangle. Note that this represents a new constraint upon mod- 
els for the structure of the emitting region surrounding Sgr A*, 
independent of those associated with the visibility amplitudes 
alone. 

1 Due to the degeneracy amongst possible triangles, given N antennas there 
are only (AT— V)(N—2)/2 independent closure phases, and thus N— 1 additional 
phases must be supplied to produce the full compliment of phase information. 
Frequently, these are obtained via self-calibration techniques, a combination 
of non-linear algorithms in which the unknown phases are chosen such that 
the resulting image satisfies various physical constraints (see, e.g., Pearson & 
Readhead 1984, for a detailed summary of self-calibration techniques). 
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The Event Horizon Telescope (EHT) is a project underway 
that will extend current mm-VLBI arrays to shorter wave- 
lengths (0.8 mm), increased sensitivity, and greater baseline 
coverage, substantially improving the ability of mm-VLBI to 
study black holes on Schwarzschild radius scales. (Doele- 
man et al. 2009a). Anticipated and potential future stations 
sites include Chile (Atacama Pathfinder Experiment, Atacama 
Submillimeter Telescope and Atacama Large Millimeter Ar- 
ray; APEX, ASTE, and ALMA, respectively), Mexico {Large 
Millimeter Telescope; LMT), the South Pole {South Pole Tele- 
scope; SPT), and the IRAM telescopes in Spain {Pico Veleta; 
PV) and France {Plateau de Bure; PdB). Among these the 
longest baselines are u ~ 1.2 x 10 4 km, corresponding to a 
maximum angular resolution of X/2u ~ 11/xas at 230 GHz 
(1.3mm) and 7.5/ias at 345GHz (0.87mm). Expected in- 
creases in bandwidth and the phasing together of elements 
within Hawaii, Chile and CARMA (e.g., Weintroub 2008), will 
lead to substantial improvements in sensitivity. As a conse- 
quence, it will become possible in the near future to measure 
mm-VLBI closure phases on a variety of additional triangles 
with uncertainties considerably smaller than that of the Fish 
et al. (2011) result. 

Here, motivated by previous efforts to model the 1.3 mm- 
VLBI visibilities using physically motivated accretion models 
for Sgr A*, we compute the closure phases implied by the ra- 
diatively inefficient accretion flow models and 1.3 mm-VLBI 
visibility amplitude fits presented in Broderick et al. (2010). 
By doing so we address three immediate questions: 

1 . Is the new closure phase estimate consistent with the ac- 
cretion flow models we have considered in particular and 
radiatively inefficient accretion flow models generally? 

2. What is the strength of the constraint placed upon physi- 
cally motivated accretion flow models and the estimates 
of black hole spin by the measured closure phase? 

3. What are the strength of the constraints that will possible 
in the near future as the EHT develops? 

In Section 2 we briefly describe the accretion models and 
how we compute the closure phases. In Section 3 we compare 
the predicted closure phases with the measured values. In Sec- 
tion 4 we predict the closure phases for the EHT and compare 
these with the estimated uncertainties of the EHT. Finally, we 
summarize our conclusions in Section 6. 

2. COMPUTING ACCRETION FLOW CLOSURE PHASES 

2.1. Accretion Modeling 

We model Sgr A*'s accretion flow as a radiatively inefficient 
accretion flow, the details of which may be found in Broderick 
et al. (2010), and references therein, and thus we only summa- 
rize the model here. 

Sgr A* transitions from an inverted, presumably optically 
thick spectrum to an optically thin spectrum near millimeter 
wavelengths. This implies that near 1.3mm Sgr A* is only be- 
coming optically thin, and thus absorption in the surrounding 
medium is likely to be important. This transition does not oc- 
cur isotropically, happening at longer wavelengths for gas that 
is receding and at shorter wavelengths for gas that is approach- 
ing. Therefore, properly modeling the structure and relativistic 
radiative transfer is crucial to producing high fidelity images. 

For concreteness, as in Broderick et al. (2010), we follow 
Yuan et al. (2003) and employ a model in which the accretion 



flow has a Keplerian velocity distribution, a population of ther- 
mal electrons with density and temperature 
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and 
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respectively, and a toroidal magnetic field in approximate (/3 = 
10) equipartition with the ions (which are responsible for the 
majority of the pressure), i.e., 
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In all of these, rs = 2GM/c 2 is the Schwarzschild radius, p is 
the cylindrical radius and z is the vertical coordinate. Inside of 
the innermost-stable circular orbit (ISCO) we assume the gas is 
plunging upon ballistic trajectories. In principle the plunging 
gas can still radiate, though in practice it contributes little to the 
overall emission due to the large radial velocities it develops. 
In the case of the thermal quantities the radial structure was 
taken from Yuan et al. (2003), and the vertical structure was 
determined by assuming that the disk height is comparable to 
p. Note that all of the models we employ necessarily have the 
spin aligned with the orbital angular momentum of the accre- 
tion flow. For the regions that dominate the mm emission, this 
assumption is well justified due to disk precession and viscous 
torques, though it may be violated at large distances. 

Thermal electrons alone are incapable of reproducing the 
nearly-flat spectrum of Sgr A* below 43 GHz. Thus it is nec- 
essary to also include a nonthermal component. As with the 
thermal components, we adopt a self-similar model for a pop- 
ulation of nonthermal electrons, 



^f.nth — nth I 
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(5) 



with a power-law distribution corresponding to a spectral in- 
dex of 1 .25 and cut off below Lorentz factors of 10 2 (consistent 
with Yuan et al. 2003). The radial power-law index was chosen 
to reproduce the low frequency spectrum of Sgr A*, and is in- 
sensitive to the black hole properties due to the distant location 
of the long-wavelength emission. 

The primary emission mechanism at the wavelengths of in- 
terest is synchrotron radiation, arising from both the thermal 
and nonthermal electrons. We model the emission from the 
thermal electrons using the emissivity described in Yuan et al. 
(2003), appropriately altered to account for relativistic effects 
Broderick & Blandford (see, e.g., 2004). Since we perform po- 
larized radiative transfer via the entire complement of Stokes 
parameters, we employ the polarization fraction for thermal 
synchrotron as derived in Petrosian & McTiernan (1983). In 
doing so, we have implicitly assumed that the emission due 
to thermal electrons is isotropic, which while generally not 
the case is unlikely to change our results significantly. For 
the nonthermal electrons, we follow Jones & O'Dell (1977) 
for a power-law electron distribution, with an additional spec- 
tral break associated with the minimum electron Lorentz fac- 
tor. For both emission components the absorption coefficients 
are determined directly via Kirchhoff's law. Images are then 
produced using the fully relativistic ray-tracing and radiative 
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transfer schemes described in Broderick & Loeb (2006a,b) and 
Broderick (2006). 

Because Yuan et al. (2003) neglected relativistic effects and 
assumed spherical symmetry, it is not directly applicable here. 
For these reasons, as in Broderick et al. (2010), the coeffi- 
cients («" th , Tei n< enid were adjusted to fit the radio spectral en- 
ergy distribution (SED) of Sgr A* (see Broderick et al. 2010, 
for details on the fitting procedure). We repeated this proce- 
dure for a large number of positions in the dimensionless spin- 
inclination (a-6) parameter space 2 , producing a tabulated set 
of the coefficients («° th , T®,n® nib ) at a large number of points 
throughout the a-6 parameter space. In all cases it was possible 
to fit the SED with extraordinary precision. From the tabulated 
values, the coefficients are then obtained at arbitrary a and 8 
using high-order polynomial interpolation. 

During the 1.3mm-VLBI observations Sgr A*'s flux var- 
ied by roughly 30%. We model this as a variable accre- 
tion rate, moving the electron density normalization up and 
down. In practice, we reduced the electron density normal- 
ization by an amount sufficient to produce a total flux of 
2.5 Jy, and then multiplied the resulting images by a correc- 
tion faction during the mm-VLBI data analysis. Because the 
source is not uniformly optically thin, this is not strictly cor- 
rect, though this makes a small change to the images them- 
selves. We produced 9090 images, with flux normalized as 
described above, at a G {0, 0.01 , 0.02, 0.98, 0.99, 0.998} for 
each 6 £ {1°, 2°, ...,89°, 90°}. We then produce models with 
arbitrary position angles, £, by rotating the images on the sky. 
For this purpose we define £ as the position angle (east of north) 
of the projected spin vector. 

2.2. Computing the Closure Phases 

The intrinsic complex visibilities are obtained in the standard 
fashion: 



(6) 



where I(a,f3) is the intensity at a given set of angular coor- 
dinates. These are subsequently modified to account for the 
observed interstellar electron scattering (see, e.g., Bower et al. 
2006; Broderick et al. 2010) 3 . The closure phase for a given 
triplet of observatories is then given by 

®ijk = arg [V(Uij,Vij)] +arg [V(u jk ,v jk )] +arg[V(u kh v ki )] . (7) 

We compute the <&y£ for all spin magnitude, inclination and 
position angles considered in Broderick et al. (2010). From 
these we determine the probability density of a given $yj: 

/>(*) = J d 3 a p(a) S [^(a) - *] , (8) 

in which p(a) is the probability of a given vector spin. For our 
purposes here we consider two forms of p(a): an isotropically 
distributed spin direction with a flat prior upon the spin magni- 
tude, the same prior as that adopted in Broderick et al. (2010), 

2 The black hole angular momentum is given by aGM 2 /c, with a = 1 corre- 
sponding to a maximally rotating Kerr spacetime. We define 8 such that models 
with = 0° are viewed along the spin axis, and thus the thick accretion disk 
viewed face-on. Conversely, models with 8 = 90° are viewed perpendicular to 
the spin axis, and thus the thick accretion disk is viewed edge-on. 

3 Note that the empirically measured Gaussian scattering kernel does not 
produce phase shifts in the visibilities, and thus we may neglect the scattering 
in computing the closure phases. Nevertheless, we must keep this effect for the 
computation of the expected closure phase uncertainties. 
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Figure 1. SMT-CARMA-JCMT closure phases, determined at 12:20 UTC on 
Day 96 of 2009, for the accretion models considered by Broderick et al. (2010) 
as a function of dimensionless spin magnitude (a) and inclination (8) for a 
number of position angles (§). Black contour lines show the closure phase 
limit of 0° ±40° reported by Fish et al. (2011). For reference, the posterior 
probability contours (white) associated with the 1-<t (solid), 2-cr (dashed), and 
3-cr (dotted) cumulative probabilities as determined in Broderick et al. (2010) 
are overlaid. Note that these are degenerate to $ — > — <E>. 

and the posterior probability of a given spin, defined in Broder- 
ick et al. (2010), after fitting the 1 .3mm-VLBI amplitudes. We 
will refer to these as the "isotropic" and "amplitude-fitted" pri- 
ors, respectively. Interpolation errors associated with the finite 
density of points in the spin parameter space and the finite size 
of the imaged region, resulted in some small-amplitude saw- 
tooth features in the As a consequence, we smoothed 
the /?($) shown on scales of l°-9°, depending upon the par- 
ticular case, with 5° being typical. Where large spike features 
appeared we did no smoothing, and in all cases the wings of 
the probability distribution remained unchanged. 

Note that fitting the visibility amplitudes necessarily imparts 
a 180° degeneracy in the position angle, corresponding to a 
sign change in the resulting closure phase. As a result, in prin- 
ciple, all closure phase estimates we present are unique only up 
to a sign. Nevertheless, we show the closure phases associated 
with one particular choice for the position angle, specifically 
that for which the spin of the most likely configuration is ori- 
ented -52° east of north. 

2.3. Closure Phase Uncertainties 

The precision with which a given closure phase may 
be measured depends upon the signal-to-noise ratio (SNR) 
of the visibilities on the individual baselines, sa = 



Vij ^2Bt /SEFD/SEFDj, where V u is the visibility amplitude, 
B is the bandwidth, r is the atmospheric coherence time, and 
SEFD, is the system equivalent flux density for a given sta- 
tion. Here we take r = 10 s, typical of the atmospheric coher- 
ence times at the wavelengths and sites of interest. However, 
because closure phases are insensitive to atmospheric phase 
fluctuations, they may be coherently averaged over time scales 
much larger than r, with the measurement uncertainties de- 
creasing as JrjT . This is limited by the time scales on which 
the baseline orientations and the intrinsic structure of Sgr A* 
change significantly. In quiescence, the former is the limit- 
ing factor, and thus we choose T = lOmin. However, during 
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Figure 2. Probability distribution of the closure phases associated with the ac- 
cretion flow models considered by Broderick et al. (2010) for the SMT-JCMT- 
CARMA triangle during the times at which Fish et al. (201 1) report measure- 
ments of this quantity. The dark grey and light grey regions are the l-cr and 
2-ct regions, defined by cumulative probability, and the light grey hatched re- 
gion shows the measured value (0° ±40°). The vertical red line indicates the 
closure phase associated with the most likely accretion flow model. For ref- 
erence, the distribution of closure phases when the 1.3mm-VLBI prior is not 
applied is shown in the upper inset. Note that these are degenerate to <I> — > — #, 



Table 1 

Estimated System Equivalent Flux Densities at 230 GHz 



Facility 



Diameter" (m) 



SEFD b (Jy) 



Hawaii 

CARMA 

SMT 

LMT 

APEX 

ALMA 



1 
1 
1 

10 



23 
27 
10 
32 
12 
38 



4900 
6500 c 
11,900 
10,000 d 
6500 
500 



a Effective aperture when number of antennas (N) are phased together. 
b Expected system equivalent flux density values toward Sgr A* in- 
clude typical weather conditions and opacities. 
c Based upon recent observations. 

d Completion of the dish and upgrades to the surface accuracy and 
receiver will eventually lower the SEFD by more than a factor of 10. 

periods of flaring activity, or if Sgr A* exhibits some dynam- 
ical structure, this could be considerably shorter. Finally, un- 
less otherwise stated, we assume B = 4 GHz, the anticipated 
near-term bandwidth target of the EHT, and adopt the SEFD, 
reported in Table 1 of Doeleman et al. (2009b), partially re- 
produced in Table 1 here. 4 Specifically we assume that the 
SMA, CSO and JCMT are phased together in Hawaii, six of 
the 10.4m and two of the 6 m CARMA antennas are phased to- 

4 This choice of bandwidth and set of SEFD, results in baseline SNRs 
approximately an order of magnitude larger than those obtained during the 
1.3mm- VLBI observations reported in Fish et al. (201 1). 



gether, and ten of the 12 m ALMA antennas are phased together 
in Chile. In terms of the Sy, the closure phase uncertainty is 
given by, 



where 



L{s). 



r 1 [Usij)Us jk )UsH)] , 



74 



(9) 



(10) 



in which Iq and I\ are the hyperbolic Bessel functions of the 
first kind (Rogers et al. 1984, 1995). In the high SNR limit 

(Sij, Sjk, s^i 3> 1) this reduces to 



r -2 -2 

j V s ij + S jk + S k 



-1/2 
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though we will make use of the exact expression. Note that 
since cr$ depends upon the individual SNR values, it is model 
dependent; the values we quote are associated with the most 
likely accretion model, as determined via the 1.3mm-VLBI 
visibility amplitudes alone. 

3. CLOSURE PHASE ON THE SMT-CARMA-JCMT TRIANGLE: 
COMPARISON WITH MEASURED VALUES 

There are two 15 min periods over which Sgr A*'s closure 
phase was measured on the SMT-JCMT-CARMA triangle on 
the night of April 6, 2009 (corresponding to day 96), beginning 
on 12:20 UTC and 12:40 UTC (Fish et al. 2011). The distri- 
bution of the closure phases at 12:20 UTC in spin magnitude 
(a), spin inclination (9), and spin position angle (£) are shown 
in Figure 1, with the closure phases at 12:40 UTC differing by 
only a small amount. For reference the l-cr, 2-er, and 3-cr pos- 
terior probability contours from the 1 .3 mm- VLBI visibility fits 
obtained in Broderick et al. (2010) are overlaid, shown by the 
solid, dashed, and dotted white contours, respectively. 

It is immediately apparent that small closure phases are not 
generic. Assuming the isotropic prior, fewer than 49% of the 
accretion flow models are consistent with the closure phase 
measurement. The inset in Figure 2 bears this out; while small 
closure phases are marginally preferred, a large closure phase 
island exists with a substantial likelihood of any phase being 
measured. The excluded models are not confined to a particular 
region in the spin parameter space, extending over a wide range 
of spin magnitudes, inclinations and position angles. Thus, de- 
spite the large uncertainties, the measured closure phase has 
considerable diagnostic ability. 

However, this changes dramatically when we restrict our at- 
tention to those models which are consistent with the 1.3mm- 
VLBI amplitudes. The white contours in Figure 1 show the 
posterior probability of a given spin vector based upon the 
1.3mm- VLBI visibility amplitudes alone, i.e., the amplitude- 
fitted prior. The regions of high probability are confined to 
a narrow valley of low closure phases. Hence, while accre- 
tion models generally are not necessarily consistent with the 
observed closure phases, the vast majority of those that were 
found to reproduce the visibility amplitudes are broadly con- 
sistent with the recent measurement. Once the amplitude-fitted 
prior is adopted, the probability of finding a closure phase that 
is larger than 40° falls to less than 3%. 

The amplitude-fitted models that violate the recent closure 
phase measurement are confined to the 2-er significance re- 
gions, primarily at large spin magnitude and large negative po- 
sition angles, where the closure phase is ~ 180°. This behavior 
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Figure 3. Probability distributions for the closure phases at 12:20 UTC on day 96 of 2009 associated with the accretion flow models considered by Broderick et al. 
(2010) for all triangles constructed from the SMT, Hawaii, CARMA, LMT and Chile stations. The hatched region in the upper left indicates the closure phase 
measurement reported in Fish et al. (201 1). In all plots the dark and light grey regions are the 1-cr and 2-cr regions, denned by cumulative probability, the vertical 
red lines show the closure phase of the most likely accretion flow model, and the red horizontal error bar indicates the accuracy with which the EHT should be able 
to measure the closure phases in the future on 10 min time scales (using 10 ALMA antennas in Chile, employing the recently phased sites in Hawaii, and 8 CARMA 
antennas). For reference the u—v tracks are shown in the upper right of each plot, with the triangle for the closure phase plotted indicated explicitly. Note that these 
are degenerate to <E> — > — <&, 
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is similar to that exhibited by annular models, in which the clo- 
sure phase varies between 0° and 180°, depending upon the an- 
gular size of the annulus. As a result, improvements of factors 
of a few in sensitivity will allow closure phase measurements 
that can compete with the visibility amplitudes in their power 
to constrain the structure of Sgr A*. 

These conclusions are borne out by the closure phase prob- 
ability distributions shown in Figure 2, computed with the 
amplitude-fitted prior. For these we show both times, which 
exhibit slightly different distributions, though in both cases the 
closure phases are dominated by peaks near -10° and 10°. 
The probability of measuring a closure phase within the al- 
lowed range exceeds 97% at both times. The closure phase 
of the model with the highest amplitude-fitted prior probabil- 
ity (based upon the 1.3mm-VLBI visibility amplitude) is lo- 
cated at 12.6°. However, significant probabilities extend out 
to <I> ~ ±30°, and thus factor of few improvements in sensi- 
tivity necessarily result in a significant likelihood of detecting 
non-trivial closure phases. 

4. CLOSURE PHASES ON FUTURE TRIANGLES 

The development of the EHT will be characterized by the 
inclusion of additional baselines and improvements in station 
sensitivity. Thus, we also compute the closure phases along tri- 
angles including intermediate and large north-south baselines. 
Figure 3 show the probability distribution for closure phases on 
all triangles made with the SMT, CARMA, Hawaii, the LMT, 
and Chile, as seen at 12:20 UTC on day 96 of 2009. Closure 
phase distributions on other days and times are qualitatively 
similar. While the closure phase on the SMT-CARMA-Hawaii 
triangle is small, this is not generic. Closure phases with mag- 
nitudes of ^45° are likely on five of ten triangles, including 
the SMT-Hawaii-Chile and Hawaii-Chile-LMT triangles, for 
which the most probable model has closure phases comparable 
to -90° (shown by the red vertical lines in each panel) 5 . 

More importantly, with improvements in the station sensi- 
tivity, arising from increases in bandwidth and the phasing to- 
gether of multiple elements at individual stations, substantial 
improvements in the precision with which closure phases can 
be measured are possible. Despite the small visibilities associ- 
ated with the most probable model (see Broderick et al. 2010), 
on all baselines constructed from the five stations we consider, 
Sgr A* is detectable on the typical atmospheric coherence time 
(10 s) with SNR exceeding unity. Over the dynamical time 
scale of Sgr A*, roughly 10 min, the uncertainty in the closure 
phase is reduced by an additional order of magnitude, resulting 
in the red, horizontal error bars in the individual panels of Fig- 
ure 3, ranging from 0.6°-l .8°. In all cases this is much smaller 
than the presently allowed range of values, and provides suffi- 
cient precision to distinguish likely values from zero. 

5. IMPROVING THE CONSTRAINTS UPON SGR A*'S SPIN 

The high precision with which closure phases may be mea- 
sured in the near future, warrants a brief examination of their 
ability to help constrain the spin of Sgr A*. To estimate this we 
consider the effects of including a single closure phase mea- 
surement on each of the various station-triangles discussed in 

5 Note that in 30% of the cases shown the red line lies outside the l-cr re- 
gion. This is not unexpected given the definition of the l-cr and 2-<r regions, 
and does not indicate that the model with the highest posterior probability is 
inherently unlikely. These deviations are generally largest for large open tri- 
angles, associated with the intrinsically larger value and distribution of closure 
phases in these cases. 



Section 4. In practice, where many closure phase measure- 
ments for each triangle are anticipated, these are expected to 
be even more constraining. 

Figure 4a shows the posterior probability distribution ob- 
tained from analysing the 1.3mm-VLBI amplitudes alone, as 
a function of spin magnitude (a) and inclination (8), after 
marginalized over the spin position angle (£) 6 . Including the 
closure phase measurement reported in Fish et al. (201 1) elim- 
inates some of the high-spin, high-inclination models, produc- 
ing the probability distribution seen in Figure 4b. Nevertheless, 
the l-cr and 2-cr contours are not significantly affected. The 
represents the present constraints upon the spin of Sgr A*. 

To illustrate the improvement associated with higher- 
precision measurements and additional triangles, we will as- 
sume the most likely model, as defined by the 1.3mm-VLBI 
visibility amplitudes alone, is the "true" model for Sgr A*'s 
emission region. This allows us to compute the visibilities 
throughout the u-v plane, obtain the SNRs at each point, and 
thus infer both the expected closure phase, $, ;J t,o and the pre- 
cision with which it may be measured, cr$, /t0 . Based upon this, 
we can construct a likelihood associated with these observa- 
tions, 



£(a) = C£| V |(a)exp 



■E 



ijk.0- 



2oi 



(12) 



where C is some constant and C\y\ is the likelihood of observ- 
ing the measured visibility amplitudes alone. This may then 
be used to produce a posterior probability density, p{a)), as de- 
scribed in Broderick et al. (2010), the particulars of which de- 
pend upon which closure phase measurements we choose to in- 
clude and their associated measurement uncertainties. In what 
follows we consider a handful of possible measurements that 
are relevant for upcoming observations. 

Improving the precision of the closure phase measurement 
on the SMT-CARMA-JCMT triangle to ±5° substantially im- 
proves the spin estimation in two respects. Firstly, it would 
conclusively eliminate the 180° degeneracy in the position an- 
gle resulting from fitting visibility amplitudes alone. Secondly, 
as seen in Figure 4c, it begins to significantly modify the l-cr 
and 2-cr contours. Thus, even with modest precisions, closure 
phase measurements can substantial improve black hole spin 
estimation in Sgr A*. 

For comparison with ongoing observations, we show in Fig- 
ure 4d the constraints associated with closure phases measured 
with a phased Hawaii, phased CARMA, SMT, and APEX, all 
with a 1 GHz bandwidth (in contrast to the 4 GHz bandwidth as- 
sumed elsewhere). Even without the substantial sensitivity im- 
provements associated with large bandwidths and phasing up 
multiple ALMA antennas, the resulting constraint is strength- 
ened dramatically, in principle allowing a spin magnitude esti- 
mate with precision ±0.1. 

The spin constraints associated with observations with a 
phased ALMA are shown in Figure 4e, at which point the al- 
lowed model parameter space becomes too small to resolve 
with the parameter sampling we employed. When the LMT 
is included (Figure 4f), the allowed region is reduced to a sin- 
gle grid point in our model space (Aa = 0.01, A9 = A£ = 1°). 
Whether such high precision spin estimation is possible in 
practice is almost certainly going to be limited by the system- 

6 Note that the allowed region is somewhat broader than shown in Figure 1 
due to this marginalization. 
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Figure 4. Marginalized posterior probability distributions for the spin magnitude and inclination, and the spin position angle (inset), for various closure phase 
measurements. In particular, shown are the marginalized posterior probabilities when: (a) only the visibility amplitudes are considered (note that this is somewhat 
more broadly distributed as a result of marginalizing over £), (b) the closure phase measurement reported in Fish et al. (201 1) is included, (c) a single closure phase 
measurement of 12.6° ±5° on the SMT-CARMA-JCMT triangle is included, (d) closure phases are measured on all triangles connecting the Hawaii, CARMA, 
SMT, and APEX stations, comparable to a single ALMA dish, assuming the presently most likely model, (e) closure phases are measured on all triangles connecting 
the Hawaii, CARMA, SMT, and Chile stations, using 10 ALMA dishes, assuming the presently most likely model, and (f) when all triangles connecting the Hawaii, 
CARMA, SMT, Chile and LMT stations are considered, assuming the most likely model. In the final instance, the relevant probability region is no longer resolved, 
dominated overwhelmingly by a single pixel. See the text for more details. 
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atic uncertainties in the modeling of the emission region. Nev- 
ertheless, the constraining power of even a handful of closure 
phase measurements is clear. 

6. CONCLUSIONS 

The measurement of mm-VLBI closure phases in Sgr A* 
provide an independent, significant constraint upon the struc- 
ture of the emission region. Models based upon radiatively 
inefficient accretion flows do not generally have small closure 
phases, and thus are not generally consistent with the recently 
measured value on the SMT-CARMA-JCMT triangle. As a 
consequence, even in the presence of large uncertainties, clo- 
sure phases can place significant constraints upon the structure 
of Sgr A*. 

However, when only accretion flow models that fit the 
1.3mm-VLBI visibility amplitudes are considered, the proba- 
bility of finding an excluded closure phase falls to less than 3%. 
Hence, the vast majority of the acceptable models in Broderick 
et al. (2010) are broadly consistent with all current 1.3mm- 
VLBI constraints. This is surprising since the 1.3 mm-VLBI 
amplitudes and closure phase measurements provide indepen- 
dent constraints, and thus having fit the accretion models to 
the former there is no a priori reason to expect that they would 
be consistent with the latter. Based upon the models that are 
consistent with the 1 .3mm-VLBI visibility amplitudes, we pre- 
dict that the closure phase on the SMT-CARMA-JCMT trian- 
gle during the time it was measured was between ±30°, with 
±13° most preferred. Improving the SNR of mm-VLBI obser- 
vations by a factor of a few, achievable by phasing up multiple 
CARMA antennas and employing the recently phased sites in 
Hawaii, should allow detection of a non-zero closure phase on 
the SMT-CARMA-JCMT triangle. Note that this would elimi- 
nate the existing 180° ambiguity in the orientation of Sgr A*. 

If Sgr A* is well described by radiatively inefficient ac- 
cretion flow models, currently acceptable closure phases on 
large triangles (e.g., SMT-Hawaii-Chile and Hawaii-Chile- 
LMT) can be as large as 90°, and typical values on triangles 
which include long baselines of 45°. Future improvements in 
sensitivity associated with the near-term development of the 
EHT will substantially increase the precision with which clo- 
sure phases can be measured. On all triangles we considered 
(constructed from the SMT, CARMA, Hawaii, the LMT, and 
Chile) the resulting precision should be sufficient to measure 
the deviations from 0°. 

Improved closure phase measurements and baseline cover- 
age result in dramatically improved constraints upon the spin 



of Sgr A*. Including a single antenna in Chile (e.g., APEX or 
a single ALMA dish) produces spin estimates with precisions 
of roughly ±0.1. With a partially phased ALMA station, this 
improves to ±0.03, and with the LMT is better than ±0.01, 
at which point we fail to resolve the likely region. In the lat- 
ter cases, the spin estimation will almost certainly be limited 
by systematic uncertainties in the modeling of the emission re- 
gion motivating the consideration of more sophisticated accre- 
tion and outflow models. 
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